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Abstract 



Tunneling through time-dependent potentials is of relevance to a number of physical 



O . 

Q , problems. Using a WKB analysis Azbel' has recently studied the effects of a nonharmonic 



time-dependent perturbation embedded in an opaque potential barrier. He suggests the 
c5 ' existence of three different regimes for transmission in such systems: direct tunneling, acti- 

vation assisted tunneling, and elevator resonant activation. We address the same problem 
Q I with a numerical technique based on wave-packet simulations. Our numerical results are 

in qualitative agreement with Azbel's picture. The outcome depends on the characteris- 
tic time T of the nonharmonic potential. There is a transition from direct tunneling to 
the activation regime around a "crossover" value Tc which is determined by the Biittiker- 
Landauer time, the distance between the initial energy and the lowest resonant level, and 
the amplitude of the time-dependent perturbation. The total transmission probability is 
strongly enhanced when entering the activation regime. 



PACS. 73.40.Gk 



Tunneling through time-modulated potential barriers has attracted considerable inter- 



est over the last ten years [|ll-0. As argued in Refs. ||10|, |11[] it may be important in a 
variety of physical problems, such as tunneling chemical reactions, charge exchange be- 
tween impurity centers and resonant tunneling in semiconductors. Until recently mostly 
harmonic time-dependent potentials have been studied theoretically. In that case the tun- 
neling particle with initial energy Ei may absorb or emit energy quanta corresponding to 
the modulation frequency u, and the transmitted and reflected beams develop "sidebands" 
at energy Ei±huj, Ei±2huj, .... Expressions for the sideband amplitudes have been found 
analytically for simple rectangular barriers 0, and formal generalizations to barriers of ar- 
bitrary shape were derived in Ref. 0]. In addition, numerical calculations involving wave 
packets have confirmed the analytic results and completed our understanding of the physics 
involved @,0, 1,1. 



In two recent papers |]10|, |T^ Azbel' showed that tunneling and activation in a nonhar- 
monic time-dependent potential can reveal new and interesting physics. Azbel' considered 
an opaque static barrier V{x) with an adiabatic time-dependent component Vrit) localized 



to a narrow region in the interior of the static part |jT2|], see Fig. 1. The potential Vrit) 



varies on a time scale T which limits the range of activation energies to roughly h/T. If 
the potential varies very slowly, there is virtually no activation, and the transmission prob- 
ability is essentially unchanged from its stationary value. However, by decreasing T there 
may, at a given instant of time to, be sufficient energy available to activate the particle to 
the lowest instantaneous resonant level Erito) in the well created by the time- dependent 
potential. Since the resonant-state lifetime r^ is long for opaque barriers, the activated par- 
ticle will be trapped in the resonant level and only escape after it has been lifted above the 
static-barrier edge by the time-dependent potential. The resonant level acts like an elevator 
for the tunneling particle; hence the name "elevator resonant activation" (ERA) [^. Ac- 



cording to Azbel' the transition from ordinary stationary tunneling to ERA happens when 
T ~ tbl, where tbl is the Biittiker-Landauer time [|ll for tunnehng out of the well with 



energy Ei. With an opaque barrier, as in Fig. 1, one has tbl — Jo° dxJm* /2[V{x) — Ei\. 
Here V{x) is the static potential barrier, xq is the classical turning point of the potential, 
i.e., V{xo) = Ei, and m* is the effective mass of the tunneling particle. For a rectangular 



barrier with V{x) = Vq and Xq = d one has tbl — dJm* /2{Vo — Ei). Azbel' also predicts 
the existence of an intermediate regime, "activation assisted tunneling" (AAT), where the 
particle is activated to an energy E < -Er(^o) before it tunnels out. Since all transmission 
processes require tunneling into the well at energy E^, the transmission probability is al- 
ways exponentially small. However, the total transmission probability is predicted to be 
greatly enhanced in the transition from direct tunneling to AAT and ERA. 

In the present work we carry out numerical simulations for tunneling systems with 
parameter values appropriate to semiconductor heterostructures. Such systems, with 
their large flexibility and tunability, seem to be a promising candidate for observing the 
new effects predicted by Azbel'. We study the tunneling dynamics by solving the time- 
dependent Schrodinger equation numerically [|1^, and we restrict ourselves to one spatial 
dimension. The initial state ^(x; t = 0) is a quasi-monoenergetic minimum uncertainty 
Gaussian wave packet with mean energy Ei, energy width AE, and mean velocity fj. 
Since we have semiconductor heterostructures in mind, suitable units of measure will be 
Angstrom [A] for length, millielectronvolts [meV] for energy, and femtoseconds [fs] for 
time. To be specific, we shall assume GaAs contacts and AlGaAs barriers with height 
230 meV. We choose the spatial width of the wave packet to be Ax = 1000 A. Then 



AE ~ {dE/dk)Ak = J2Ei/m*h/Ax ^ 5 meV if Ei equals half the barrier height, and 
m* = 0.067mo (the effective mass in GaAs). Thus the initial wave packet has a well-defined 
energy in all the numerical examples below, in the sense that AE{dPtr/dE)/Ptr << 1 P 



Here Ptr is the instantaneous transmission probability of the potential barrier at energy 
Ei. The initial state has its center of mass < a;(0) > more than 4000 A away from the 
barrier structure. Thus the overlap between il){x; t = 0) and the barriers is negligible. 

Our choice of initial state is not of purely academic interest. By using modulation- 
doped semiconductor heterostructures it is possible to fabricate tunneling barriers which 
transmit a fairly monoenergetic beam of electrons. These electrons can move ballistically 
over tens of nanometers [ffH] and should be well described by wave packets of the kind used 



here. 

The static part of the potential V{x) can be made by standard epitaxial-growth tech- 
niques. Alternatively, the transport could take place in a two-dimensional electron gas. In 
that case the potential-barrier structure is created by applying external gate potentials. 
Possible ways of producing the perturbation VT{t) could be by means of ultra-short laser 
pulses or time-dependent gate potentials [jl6|. It may be difficult to control Vxit) on the 



length and time scales which are typical for the systems that we have in mind (of the order 
of 100 A and 100 fs, respectively), although laser pulses with a duration less than 100 fs 
have already been demonstrated in several experiments. 

We will present numerical results for two different kinds of potential-barrier structures, 
see Fig. 2. Structure A is a single rectangular barrier of width Wq and height Vq. The 
time-dependent part of the potential acts on a segment of width Wi located at the center 
of the static barrier. In this region the total potential is Vq when \t\ ^>- oo and Vq — Vi 
for t = to. We have used Vi = Vq in all the examples below. Within the adiabatic picture 
the lowest instantaneous resonant level starts at Eoo = Er{\t\ -^ oo) ( i.e., the "virtual" 
state above the barrier), falls down to Eq = Er{to) at t = to, and rises back to E^q as 
t — > oo. Structure A will be used to investigate ERA. Structure B is similar to A, but has 
in addition a wide and low barrier of height V on each side, leading to a total width W 



for the structure. There is still a resonance in the time-dependent well, moving from Eoo 
to Eq and back again. The positions of Eq and E^o can be tuned by changing the well 
width Wi (cf. Figs. 2B1 and 2B2). In addition there are one or more "quasi- resonant" 
levels Ej > V {j = 1,2, ...) corresponding to resonances above a single rectangular barrier 
of width W and height V. The positions of these levels are quite insensitive to both the 
instantaneous value and the width of the time-dependent potential. With structure B we 
will demonstrate AAT and ERA, where the outcome will depend on the relative positions 
of Eq and the various Ej. In order to observe a transition from direct tunneling to the 
activation regime, the initial-state energy Ei must be chosen to lie below all the resonant 
levels. In the opposite case, with Eq < Ei, the elevator effect will be present for any finite 
T since no activation energy is required to trap the particle in the resonance |TI|]. In the 



tables in Fig. 2 we have collected numerical values for the most important parameters in 
the examples described below. We have also plotted the transmission probability of the 
instantaneous potential, both at t = to and t ^ oo. Notice the peaks in Ptr for E = Ej in 
Fig. 2B. Although the amplitudes of these peaks depend strongly on the value of Vxit), 
their positions remain more or less unchanged. The parameter to is chosen such that the 
center of mass at t = tg, < a;(to) >=< x{0) > +fjto, of a freely moving wave packet with 
group velocity Vi would coincide with the center position of the time-dependent part of 
the potential. In this way our simulations are as close as possible to the stationary case 



(Ax -^ oo) studied in Refs. ||T0|, [TT 



The exact form of the time-dependent part of the potential is not important. However, 
in order to observe ERA, it is crucial to have a nonharmonic variation in time. A nonhar- 
monic time dependence implies, as we shall see below, a continuous activation probability. 
Under these conditions a particle which has been activated from E^ to the lowest instan- 
taneous resonant level Eq, can continuously absorb infinitesimal amounts of energy and 



follow the resonance as it rises from Eq to its maximum value Eoq. In contrast, in a purely 
harmonic potential with frequency u a particle can only absorb energy in multiples of hu, 
and hence cannot be trapped in the continuously moving resonant level. In the numerical 
calculations below we will use the time-dependent potential 

Vi 
^^^'^ = "cosh((t-io)/T)' ^^^ 



which was also used in an example in Ref. |10| 



Let us try to gain some insight into when the transition from direct tunneling to acti- 
vation takes place. This must happen when the characteristic time T of the nonharmonic 
potential is such that the probability Pact{Ei — > Eq]T) for activation from Ei to the lowest 
resonant level Eq is approximately equal to the probability Pdir{Ei) for direct tunneling 
out of the time-dependent region with the initial energy E^. An estimate of the activation 
probability can be made by using standard time-dependent perturbation theory [0. We 
find 



h J—oo 

ttViT 



(2) 



hcosh[l7TT{Eo-Ei)/h] 

where \i > and |0 > denote the initial state and the lowest resonant state, respectively, 
and we have used Vr(t) from Eq. (1). We have also taken into account that Pact is the 
activation probability conditional on the particle having tunneled into the time- dependent 
region, and we assume that the lowest resonant state can be approximated by a symmetric 
and normalized wave function. A "crossover" time Tc can now be defined via 

Pact{Ei^Eo-T,)c^Pdir{Ei). (3) 

Since the barrier is opaque when seen from the time-dependent region, we have 



Pdir{Ei) ~ exp[— 2(iw2m*(\/o — Ei)/h], and hence 



ixV.T, ^' 



eM-2dJ2m*{Vo-E{)/h]. (4) 



hcosh.[\'KT^{EQ-Ei)/?i]_ 
We have found that the argument of the cosh-function in Eq. (4) is typically much larger 
than unity. In that case we may rearrange Eq. (4) and obtain a simpler equation for the 
crossover time, 

■K Eq- Ei 7r(Eo - Ei) h 

From Eq. (5) it is evident that the Biittiker-Landauer time tbl sets a scale for the crossover 

from direct tunneling to activated processes. As mentioned earlier, this was pointed out 



very clearly in Ref. |10|. However, equally important parameters in this connection are the 
distance from the initial energy to the lowest resonance, Eq — E^, and the amplitude of the 
potential modulation, Vi. 

In order to analyze the tunneling and activation process in energy space, we perform a 
spatial Fourier transform of the wave packet in the asymptotic limit t ^ oo, thus obtaining 
the momentum (or energy) distribution of the transmitted {p > 0) and reflected {p < 0) 
parts of the wave packet [|18 . 



In our first example we use the structure from Fig. 2A. In Fig. 3 we have plotted 
on a logarithmic scale the energy distribution of both the transmitted (tr) and reflected 
(ref) parts of the wave packet for a fairly slow potential variation (T = 182 fs). One can 
see that most of the wave packet is reflected at energy Ei without ever reaching the time- 
dependent region. In other words, we are in the opaque barrier regime, consistent with the 



assumptions in Refs. [0, |Tl|]. Next, we see that for this slowly varying potential the energy 
distribution of the transmitted part of the wave packet is also centered at Ei, indicating 
that hardly any activation has taken place. Notice, however, the small peak (or "shoulder") 
at Eoo, i.e., at the position of the resonant level as t — ;► oo. These energy components are 



activated to Eq by frequencies u ~ {Eo — Ei)/h » 1/T in the exponential tail of the 
Fourier spectrum of Vrit) and subsequently follow the resonance adiabatically from Eq to 
Eoo- Finally, notice the symmetry of activated energy components in the reflected and 
transmitted parts of the wave packet. This is seen in all our numerical examples because 
activated components of the wave packet "see" a spatially symmetric potential. Hence 
they have equal probability of tunneling out in either direction and thus contribute with 
the same amount to reflection and transmission. 

In the following we present numerical results which illustrate the transition from direct 
tunneling to activation. We will focus on the transmitted part of the wave packet, and in 
Fig. 4 we have plotted its energy distribution (now on a linear scale) for selected values of 
T, and for the three different structures presented in Fig. 2. 

The first panel in Fig. 4A is a linear plot of the transmitted distribution in Fig. 3, and 
it is clear that direct tunneling dominates completely for such a slowly varying potential. 
Thus the total transmission probability for the wave packet PwpiT = 182 fs) is very little 
different from the stationary value Pwp{T — > oo) (see the rightmost panel in Fig. 4A; solid 
dot and dotted line, respectively) [|1^. A faster potential, T = 22 fs, yields a strikingly 
different picture. The energy spectrum of the transmitted packet is now dominated by 
components around E^q. In this case ERA dominates, and there is indeed a huge increase 
in the transmission probability. Finally, we have included the case with T = 91 fs, which 
represents an intermediate situation with roughly equal amounts of direct tunneling and 
ERA. This is in good agreement with Eq. (5), which predicts a crossover between direct 
tunneling and ERA at Tc ~ 80 fs. 

As we decrease T below 20 fs, the total transmission probability goes through a max- 
imum value and then falls off rapidly. The reason is that the time-dependent potential, 
and hence the lowest resonance, is sufficiently low to permit "effective" activation only in 
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a period of time of the order of T. In other words, for decreasing T a greater part of the 
wave packet is directly reflected [^]. By varying the initial energy Ei we have further 
verified that the onset of ERA scales as expected from Eq. (5). One should also bear in 
mind that when T^IO fs, we are no longer in the adiabatic regime h/T << Ei, Vq — E^, 



and thus outside the limits of the analytical work in Refs. [T^, [TTJ- 

Activation assisted tunneling, AAT, implies activation to an energy below Eq before 
tunneling out of the time-dependent region. In our numerical calculations this should 
show up as a significant weight between Ei and Eq in the energy spectrum of the trans- 
mitted wave packet. With a structure of type A (and with an activation probability as 
above) we have seen no indication of AAT. This is not surprising since both the (static) 
transmission coefficient and the activation probability behave monotonically as functions 
of energy between Ei and Eq. Thus there is no energy between Ei and Eq at which the 
particle is preferably transmitted, and one has a transition directly from direct tunneling 



to ERA when T is decreased. We should point out that this was predicted in Ref. |jTO| for 



a similar time-dependent potential, but with an inverted parabola static part instead of 
the rectangular shape used here. 

In order to observe AAT we will in the next examples use a structure of type B (see Fig. 
2B). As discussed earlier such a potential has quasi- resonances at Ej > V to which the 
tunneling particle may be activated. Since the quasi-resonances Ej do not move with the 
time-dependent potential, the activated particle will eventually tunnel out with energy Ej, 
and one has AAT rather than ERA. Again we use a symmetric potential, with parameters 
given in the tables of Fig. 2B. The time dependence is, as above, given by Eq. (1). In 
practice the low (and wide) barriers could be made from Al^.Gai_2.As, with a lower Al 
content x than in the high (and narrow) ones. By varying the width Wi of the time- 
dependent region, the lowest time-dependent resonance Eq is tuned to different positions 



relative to the quasi-resonances Ej. In Fig. 4B1 we use a narrow well {Wi = 18 A). As 
a result Eq lies above Ei and E2, but below E-^. In Fig. 4B2 the well is somewhat wider 
(Wi = 35 A), and Eq lies below all Ej. We have again plotted the energy distribution of 
the transmitted part of the wave packet for various values of T. In the case of the narrow 
well the direct tunneling component is significant for all T. This shows up as a peak at 
E' ~ £^4 in the transmitted energy distribution. Furthermore the selected examples in Fig. 
4B1 illustrate A AT to Ei when T = 91 fs, AAT to Ei and E2 (the latter only as a weak 
shoulder in the energy distribution) when T = 46 fs, and a combination of AAT to Ei and 
E2 and ERA from Eq to E3 and above when T = 6 fs. With the wider well we have observed 
substantial activation to Eq followed by ERA already for T ~ 350 fs. When T = 182 fs 
(first panel, Fig. 4B2), most of the activated components follow the resonance to Ei and 
tunnel out. Reducing T to 46 fs results in ERA to Ei and E2. In that case the lifetime of 
the lowest quasi-resonance, fi ~ 70 fs, is so long that many of the "elevating" wave-packet 
components are trapped in the moving resonance and do not tunnel out at Ei. However, 
the second quasi-resonance has a lifetime f2 ~ 35 fs which is shorter than T. Thus the 
elevating components that pass through the level at Ei tend to tunnel out at £'2- Upon 
further reduction, to T = 6 fs, we observe ERA to Ei, E2 and E^, and some components 
even follow the resonance all the way up to E^o before they tunnel out. In Figs. 4B1 and 
4B2 (rightmost panel) we have also plotted the total transmission probability as function 
of T for the two cases discussed above. As in Fig. 4 A there is a large enhancement of P^p 



when T is made short enough to allow for AAT or ERA, as predicted in Refs. ||10| , pi]] . The 
crossover values calculated from Eq. (5) are Tc = 137 fs for Fig. 4B1 and Tc = 285 fs for 
Fig. 4B2, both in good agreement with what is observed. For sufficiently small values of 
T the transmission probability goes through a maximum and falls off rapidly, for the same 
reason as discussed in connection with Fig. 4A. 
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We have also examined other time-dependent potentials, both harmonic and nonhar- 
monic ones. As expected the qualitative features in the examples above stay unchanged as 
long as the activation probability is a continuous function of energy. On the other hand, 
with a purely harmonic potential, Vi cos cut, there is no elevator effect but simply activation 
and tunneling at energies Ei + nhuj (n = 0, ±1, ±2, . . .). 

In conclusion, we have studied tunneling through a nonharmonic time-dependent po- 
tential by letting Gaussian wave packets scatter off the barrier. Analysis of the energy 
distribution of the transmitted wave packets for suitably chosen tunneling structures con- 
firms Azbel's prediction of a crossover from direct tunneling to elevator resonant activation 
and activation assisted tunneling. An estimate shows that activation dominates over direct 
tunneling when the characteristic time of the time-dependent potential is reduced below a 
crossover value which involves the Biittiker-Landauer time, the distance from initial energy 
to the lowest resonant level, and the amplitude of the time- dependent potential. The total 
transmission probability is strongly enhanced, in our numerical calculations by one or two 
orders of magnitude, when comparing activation with direct tunneling. 

We would like to acknowledge K. Flensberg for useful discussions. 



* Also at Microelectronics Centre, Technical University of Denmark, 2800 Lyngby, Den- 
mark. 
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Figure Captions 



FIG.l. One-dimensional potential barrier V{x) with a time-dependent component Vrit) 
embedded near the center, xq denotes the classical turning point for incoming energy Ei. 
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Solid and dotted vertical lines in the well illustrate three possible processes, as described 
in the main text: direct tunneling at Ei (DT), activation to E and tunneling (AAT), 
activation to Er{to) and elevation to Er{oo) (ERA). 

FIG. 2. Potential-barrier structures, static transmission probabilities versus energy Ptr{E), 
and tables with potential-barrier parameters and resonant-level energies for the two types 
of structure used in our numerical calculations. (A) A rectangular static barrier of height 
Vq and width Wq. The time-dependent part is located at the center of the static barrier 
and has width Wi and maximum value —Vi. We have used Vi = Vq in all numerical 
calculations. The lowest instantaneous resonance moves between E^o and Eq. (B) Similar 
to the structure in A, but here we have in addition wide and low barriers of height V on 
each side. The total width is W. In addition to the moving resonance there are "quasi- 
resonances" at Ei,E2,.... Figs. Bl and B2 represent results for structures which are 
identical apart from having different well width Wi. The levels in the well labeled with 
1 and 2 denote the extremes Eq and Eoo of the moving resonance (for the narrow and 
wide well, respectively). The transmission probability is in all three cases plotted for the 
maximum and minimum value of the time-dependent potential: Solid lines are for t ^ oo 
whereas dashed lines are for t = to, i.e., for V = Vq and V = 0, respectively, in the 
time-dependent region. The peaks in Ptr are associated with the various resonant levels 

Eq, Eoo, Ei, E2, ■ ■ ■■ 

FIG. 3. Logarithm of the energy distribution |<l>(i5)p of the reflected (ref) and transmitted 
(tr) parts of the wave packet after scattering off the structure in Fig. 2A. (Note that 
the E-axis has positive values in both directions.) The characteristic time of the time- 
dependent potential is T = 182 fs. Both parts of the wave packet are dominated by energy 
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components around the initial energy i5j ^ 81 meV. The "shoulders" with their maxima 
close to -Eqo — 248 meV represent energy components that have been activated to the 
lowest instantaneous resonance Eq ~ 127 meV and elevated from Eq to E^o (ERA). 

FIG. 4. Numerical results corresponding to the structures presented in Fig. 2. Reading 
horizontally, the first three panels show the energy distribution (on a linear scale) of the 
transmitted part of the wave packet for selected values of T. Note that the scale along the 
vertical axis varies from curve to curve. The rightmost panel shows the total transmission 
probability as function of T. These curves result from integrating the transmitted energy 
distribution for each value of T. The solid dots represent the selected values of T. The 
dotted horizontal line (for the case B2 this line lies very close to P^p = 0) represents the 
stationary value of the transmission probability, Pwp{T -^ oo), for the particular structure. 
(A) Results for the rectangular barrier in Fig. 2A. The first frame (T = 182 fs) is a linear 
plot of the transmitted part in Fig. 3. (B) Results for the structures of Fig. 2B, for the 
narrow well (Bl) and the wide well (B2). In all cases the peaks in the transmitted energy 
distribution are connected to the initial energy Ei or to the various resonant levels, which 
can be identified from the peaks in the static transmission probabilities in Fig. 2. 
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